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We show that many topological and geometrical properties of complex pro- 
jective space can be understood just by looking at a suitably constructed 
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them. 
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1. Introduction. 

Most people are familiar with the sphere S^; how to make a map of a sphere 
and how to use it to get insight into the geometrical and topological properties 
of a sphere. It is rare to see complex projective space CP'^ treated in the 
same way, even though a large number of people deal with the geometrical 
and topological properties of this space in their everyday work. (It happens 
that CP*^ is precisely the space of pure states of a particle of spin n/2, so 
indeed undergraduate physics students deal with it every day! For n = 1 it 
happens that CP^ = and we get the "Bloch sphere" but for n > 1 CP" 
is not a sphere.) 

It is the purpose of this paper to convince the reader — or at least those 
readers who like this sort of thing — that we can literally draw maps of CP^ 
and CPl Now every map distorts geography, so if we choose some property 
of CP^ that we wish to illustrate it is not obvious that it will be recognizable 
on the map. Our interest lies in illustrating the geometry of quantum me- 
chanical entanglement, and the pleasant surprise is that our map works very 
well. Unfortunately we cannot really draw our map when the real dimension 
of the space to be mapped exceeds six, which means that we must restrict 
ourselves to the entanglement of two qubits in a pure state. This case is fully 
understood — it is probably the only case where this is true — so that we will 
not derive any new insights. Nevertheless we found the picture pleasing, and 
we would like to share it with others. 

The paper is organized as follows: In section 2 we introduce the picture 
and use it to illustrate some topological properties of CP^. In section 3 we 
study the separable and maximally entangled states in CP^. This is the 
lowest dimension where entanglement can occur. In section 4 we study sub- 
manifolds of states with intermediate entanglement. Section 5 is an aside 
on the Schmidt decomposition and the statistical geometry of density ma- 
trices and section 6 is another aside on the symplectic geometry of CP"'. 
Throughout we try to keep track of which properties that are special to low 
dimensions, and which are not. Our main purpose is pedagogical although 
a few new results are included (eg. in sections 4 and 6). It will be helpful 
but we hope not necessary if the reader has a nodding acquaintance with the 
3-sphere and the Hopf fibration; there are references |I| that explain all that 
is needed in elementary terms. It only remains to add that the picture was 
not invented by us. There is a branch of mathematics called toric geome- 
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try whose subject matter (roughly speaking) consists of spaces that can 
be depicted in this way. Moreover the picture has featured in the quantum 
mechanics hterature already |^ 

2. The picture. 

A pure state in quantum mechanics is described by a vector in an complex 
dimensional vector space; in Dirac's notation 

l^) = E^1«), (1) 

where \a) is a given orthonormal basis, n = N — 1 and has complex 
components. It is understood that vectors that differ only by an overall 
complex factor count as the same state. This means that there is a one-to- 
one correspondence between the set of pure states and the set of equivalence 
classes 

{Z^,Z\Z^, ...,Z'') ^ z{Z^,Z\Z^, ..^Z"") , zeC (2) 

By definition this is complex projective space CP". The numbers Z° are 
known as homogeneous coordinates. To make this description more concrete, 
suppose n = 2. Choose the complex number z and the relative phases vi and 
z/2 SO that 

(ZO,Z\z2) = (no,nie*'^Sn2e*^^) , (3) 

where < z/j < 27r and the real numbers nQ,ni,n2 are non- negative, no > 
, 77-1 > , ?T,2 > 0, and obey the constraint 

nl + nl + nl = l . (4) 

(So we work with normalized state vectors.) This means that the set of 
allowed numbers no,ni,n2 are in one-to-one correspondence with points on 
the positive octant of the 2-sphere, while the periodic coordinates z/i, z/2 are 
in one-to-one correspondence with the points on a torus. The description 
breaks down at the edges of the octant since then the torus is undefined. 
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Apart from this we already have a picture of the topology of CP^ which 
turns out to be quite useful. 

Even more to the point, this picture reflects the geometry of CP^. There 
is a natural way to define the "distance" between two pure states in quantum 
mechanics — namely in the sense of statistical distance [^. There is also a 
mathematically natural notion of distance on CP" called the Fubini-Study 
metric. The two of them coincide; the distance d between two states is given 
by 
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where Z ■ Z = Z'^Z^ and is the row vector whose entries are the complex 
conjugates of the entries of the column vector Z". Note that the maximum 
distance between two points equals 7r/2; the precise number is a convention 
but the fact that there exists a maximal distance is not. In infinitesimal form 
the Fubini-Study distance becomes the metric tensor 

,9 Z ■ ZdZ ■ dZ — Z ■ dZdZ ■ Z . . 

ds^ = = = . 6) 

z-zz-z ^ ' 

The point we are driving at is that this metric takes a very nice form in the 
coordinates given above; common algebraic work shows that 



d^ = driQ + dni + dn\ + 

(7) 

+n\{\ — n\)diJl + n\{l — n\)di'2 — 2n\n\dvidi'2 ■ 

The first piece here, given eq. (^), is recognizable as the ordinary "round" 
metric on the sphere. The second part is the metric on a flat torus, whose 
shape depends on where we are on the octant. Hence we are justifled in 
thinking of CP^ as a set of flat 2-tori "parametrized" by a round octant of a 
2-sphere. There is an evident generalization to all n. In particular for n = 1 
we obtain a one parameter family of circles (that degenerate to points at the 
end of the interval). A moment's thought will convince the reader that this 
is simply a way to describe a 2-sphere — and indeed it is well known that CP^ 
is the same thing as the 2-sphere S^, usually called the Bloch sphere when it 
is regarded as the space of states of a spin 1/2 particle. (The choice of the 
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metric means that the radius of this sphere is actually 1/2, so that the 
maximum distance between two points on the Bloch sphere is '/r/2.) 

To make the case n = 2 quite clear we make a flat map of the octant. 
Two methods suggest themselves, namely stereographic |l|] and gnomonic 
projection. The first is a standard coordinate system with several advantages; 
notably one can cover an entire sphere minus one point with one map. For 
our purposes the gnomonic (central) projection works even better; here the 
projection from the sphere to the fiat map is made from the center of the 
sphere. It does not matter that only half the sphere can be covered since 
we need to cover only one octant anyway. A decided advantage is that the 
geodesies on the sphere, i.e. its great circles, appear as straight lines on the 
map. This is obvious because a great circle is the intersection between the 
sphere and a plane through the origin, and this will appear as a straight line 
when we project from the origin. 

We choose to center the projection at the center of the octant and adjust 
the coordinate plane so that the coordinate distance between a pair of corners 
of the resulting triangle is one. Explicitly, first we choose an auxiliary basis 
in three space so that = 1 labels the center of the octant: 

X° = -^{no + ni + n2) 

X' = j={-no + n,) . (8) 

X^ = ^(-no -ni + 2^2) 

Next we define the gnomonic coordinates xi,X2 by 



Xi = ^^ ^ X' = ^^ X' = ^^- r'^xl + xl.{9) 
V6X0 yrT72 VI + 6r2 1 2 V ; 

The octant is bounded by great circles and therefore it now appears as a 
triangle centered at the origin. Its sides have coordinate length one and in 
the gnomonic coordinates the metric takes the form 



6 

ds"^ = dnl + dnl + dnl = ^((1 + Gr'^)dx ■ dx — 6(x ■ dxY) (10) 

(1 + or"')^ 

where x ■ dx = xidxi + X2dx2 and so on. Figure 1 should make all this clear. 
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Figure 1: An octant of a sphere and two flat maps: To the left we see the positive 
octant of the sphere. In the middle there is a stereographic map of the octant 
centered at (no, ni,n2) = (1, 0, 0). On the right there is a gnomonic map centered 
at -^(1, 1, 1). The projections are done respectively from the South Pole and from 
the center of the sphere, as shown. On each map four geodesies are drawn. A pair 
of these geodesies that meet an edge divides that edge into three equal parts. This 
gives an idea about how distances are distorted by the maps. 



This takes care of the octant (the "manifold with corners" in the language 
of toric geometry 0). It becomes a picture of CP^ when we remember 
that each point in the interior really represents a fiat torus, conveniently 
regarded as a parallelogram with opposite sides identified. The shape of the 
parallelogram is relevant. According to eq. (pi) the lengths of the sides are 



Li = J ds = 2'Kni\Jl — rii and L2 = 27rn2\/l — . (11) 
The angle between them is given by 

cos ^12 = , ^ = . (12) 



l-nf^l-n. 

The point is that the shape depends on where we are on the octant. So does 
the total area of the torus, 

A = LiL2sm9i2 = 47r^nonin2 . (13) 
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The "biggest" torus occurs at the center of the octant. At the boundaries the 
area of the tori is zero. This is because there the tori degenerate to circles. 
Figure 2 shows how this happens. In effect an edge of the octant is a one 
parameter family of circles, in other words it is a CP^. 

It is crucial to realize that there is nothing special going on at the edges 
and corners of the octant, whatever the impression left by the map may 
be. Like the sphere, CP" is a homogeneous space and looks the same from 
every point. To see this, note that any choice of an orthogonal basis in a 
3 dimensional Hilbert space gives rise to 3 points separated by the distance 
from each other in CP^. In projective geometry the triplet of points 
arising from an orthogonal basis in the underlying vector space is known as 
a triangle of reference. By an appropriate choice of coordinates we can make 
any such triple of points sit at the corners of an octant in a picture identical 
to the one above. 

To get used to the picture let us consider some submanifolds. CP^ is 
also known as the "complex projective plane". Every pair of points in this 
"plane" defines a unique "complex projective line" containing the pair of 
points, and such a "line" is a CP^. Conversely a pair of complex projective 
lines always intersect at a unique point. Since a CP is always a sphere (of 
radius 1/2) the terminology may boggle some minds, but these intersection 
properties are precisely what defines a line in projective geometry. Through 
every point there passes a 2-sphere's worth of complex projective lines, con- 
veniently parametrized by the way they intersect the "line at infinity", that 
is the set of points at maximal distance from the given point, which in itself 
is a CP^. This is easily illustrated provided we arrange the picture so that 
the given point sits in a corner. 

Another submanifold is the real projective plane RP^. It is defined in a 
way analogous to the definition of CP^ except that real rather than complex 
numbers are used. The points of RP^ are therefore in one-to-one correspon- 
dence with the set of lines through the origin in a three dimensional real 
vector space and also with the points of S^/Z^, that is to say the sphere 
with antipodal points identified. In its turn this is a hemisphere with an- 
tipodal points on the equator identified. RP^ is clearly a subset of CP^. It 
is illuminating to see how the octant picture is obtained, starting from the 
stereographic projection of a hemisphere (a unit disk) and folding it twice, 
as in Figure 3. 

Brody and Hughston give a beautiful account of how the physics of a 
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Figure 2: The octant really parametrizes a family of fiat tori of varying sizes and 
shapes. In a) we indicate how the torus lying over each interior point changes with 
position in the octant. The position in the octant is given by an unnormalized 
vector. At the edges the tori degenerate to circles so that the edges are complex 
projective lines, while the corners of the octant represent points. Sometimes it 
is convenient to perform a little cutting and gluing of the parallelogram before 
thinking about the shape of the torus it defines, as indicated with dashed lines 
for the torus lying over the point (1,4, 1). The size of the octant relative to that 
of the tori is exaggerated in the picture. To bring this home we show, in b), 
the largest torus — the one sitting over (1, 1, 1) — decorated with three times three 
points corresponding to three mutually unbiased bases (represented respectively by 
crosses and filled and unfilled dots). They are all mutually unbiased with respect 
to the basis that forms the corners of the octant. The coordinates {1^1, 1^2) are given 
for one basis. 
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Figure 3: Using stcrcograpliic ratlicr tlian gnomonic coordinates we show how the 
octant picture of the real submanifold RP^ is related to the standard description 
as a hemisphere with antipodal points on the equator identified. 

spin 1 particle is tied to the geometry of CP^. The space of all possible spin 
up states (with respect to some direction) forms a sphere of radius l/-\/2. 
This is not a complex projective line because of its size. The space of all 
possible spin states (with respect to some direction) forms a sphere with 
antipodal points identified since such states do not distinguish up and down, 
hence this is an RP^. This is illustrated in Figure 4 under the assumption 
that the corners of the octant correspond to eigenstates of the operator S^. 
It is interesting to observe that if we start from a CP^ (and place it as an 
edge in the picture) then we can increase the size of the sphere by deforming 
the edge, but we cannot shrink it. Therefore CP^ contains non-contractible 
2-spheres, just as a torus contains non-contractible circles. Note that two 
non-contractible 2-spheres always touch at a point, essentially because two 
complex projective lines intersect in a unique point. 

Next choose a point. Place it at a corner of the octant and surround 
it with a 3-sphere consisting of points at constant distance from the given 
point. In the picture this will appear as a curve in the octant with an entire 
torus sitting over each interior point of the curve. Readers familiar with the 
Hopf fibration know that can be thought of as a one parameter family of 
tori with circles at the ends. The 3-sphere is round if the tori have a suitable 
rectangular shape. As we let the 3-sphere in Figure 5 grow the tori get more 
and more "squashed" by the curvature of CP^, and the roundness gradually 
disappears. When the radius reaches its maximum value of 7r/2 the 3-sphere 
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Figure 4: The physics of spin systems (above). The dashed curve represents a 
sphere with radius l/\/2, and the sohd hue an RP^. Both are one parameter 
families of circles with points at one end. The sphere has a point at both ends, 
while for the real projective plane the circle sits in the torus in such a way that 
when the torus is squashed to a circle then the circle wrapped inside it suddenly 
collapses to a circle of half the size it had just before. We show how this happens 
by drawing two of the tori explicitly; again points on the sphere are represented by 
dashed lines and points on the RP^ by solid lines. For comparison we also show 
a complex projective line (below). This is a one parameter family of circles with 
points at both ends, as shown in the tori on the right hand side. Note also that 
any edge of the octant is a CP^, namely the "line at infinity" with respect to the 
point in the opposite corner. 
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Figure 5: The set of points at constant distance from a corner form a squashed 
3-sphere. In a) we show how such a submanifold appears in the octant. All the 
points in the torus lying over a point on the curve are included. In b) we show how 
the size and shape of the torus change as we move along the curve in the octant; 
at the ends of the interval the tori collapse to circles. For comparison, in c) we 
show the corresponding picture for a round 3-sphere. 

collapses to a 2-sphere, namely to the projective line "at infinity". Readers 
unfamiliar with the Hopf fibration may wish to consult Urbantke |^ at this 
point; let us just remark that only odd dimensional spheres can be squashed 
in this specific way (so it is no use to try to picture this in terms of the 
2-sphere). 

Finally, a warning: The picture distorts distances in many ways. For 
instance, the distance between two points in a given torus is shorter than it 
looks, because the shortest path between them is not the same as a straight 
line within the torus itself. Technically, the flat tori are not totally geodesic. 
Note also that we have chosen to exaggerate the size of the octant relative 
to that of the tori in the pictures. To realize how much room there is in 
the largest torus, consider mutually unbiased measurements [0. It is known 
that one can find 4 sets of orthonormal bases in the Hilbert space of a spin 1 
particle such that the absolute value of the scalar product between members 
of different basises is 1/ — clearly a kind of sphere packing problem if 
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translated into geometrical terms by means of eq. (|^). If one basis is rep- 
resented by the corners of the octant, the remaining 3 times 3 basis vectors 
are situated in the torus over the center of the octant. (See figure 2.) 

3. Separable and maximally entangled states in CP'^. 

We now wish to illustrate entanglement. This forces us to increase the dimen- 
sion, since our system should be composed of two subsystems. If we choose to 
study pairs of entangled qubits the complex Hilbert space is ® = C^, 
and the space of pure states becomes the six real dimensional space CP^. 
This is complex projective space. A preliminary remark about it is that it 
contains complex projective planes and lines (CP^ and CP^) as submani- 
folds, and the intersection properties of these are just those of planes and 
lines in ordinary Euclidean space, with the important simplification that all 
troublesome exceptional cases (such as parallel planes that do not intersect 
in a line) are absent. 

It is easy to make a picture of this six dimensional space. Choose the 
coordinates 

{Z^,Z\Z\Z'') = {n^,n,d^\n^d^\n^d^-^) . (14) 

The phases now form a 3-torus (that we can picture as a rhomboid), while 
the non-negative real numbers no etc. obey 

nl + n\ + nl + nl = l . (15) 

Hence they form a hyperoctant of the 3-sphere. Without further ado it is 
clear that if we perform a gnomonic projection centered at 

(no,ni,ri2,?^3) = ^(1, 1, 1, 1) (16) 

then we obtain a picture of this hyperoctant clS cl TG gular tetrahedron. In 
this picture straight lines correspond to geodesies on the round hyperoctant. 
Over each point in the interior there is a flat 3-torus of a definite shape. The 
faces of the tetrahedron are complex projective planes, its edges are complex 
projective lines, and its corners are points. Such pictures will appear soon. 
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We do not give the transformation to gnomonic coordinates here. This 
is because their only advantage is to give a nice symmetrical picture that is 
easy to draw, and most of our pictures can be drawn using only the fact that 
geodesies are straight lines, plus explicit knowledge of the two dimensional 
case. (Calculations should be performed in coordinate systems suited to 
calculation — for most purposes the embedding coordinates uq, rii, n2, will 
do.) 

Now, what about entanglement? Our illustrations in this section will 
depict submanifolds of states of constant entanglement, namely separable 
states (that are not entangled) and maximally entangled states. In the next 
section we show how states of "intermediate entanglement" sit in CP^. The 
same story has been told before but not in quite this way. 

A brief review of the facts seems appropriate. State vectors for composite 
systems are conveniently written as 

-| n n 

l^) = ^EEcd^)U)^ (17) 

where Cij is an x iV matrix with complex entries and during the review 
we keep the state vector normalized. Throughout the discussion we rely on a 
fixed way of splitting the Hilbert space into a tensor product of two smaller 
Hilbert spaces. In other words it has been agreed that the Hilbert space 
is TiA ®'Hb in a specific way. Otherwise the term "entanglement" has no 
meaning. For the 2x2 case let us agree that 

{Z^ , , Z'^ , Z"^) = (Coo, C*oi) C'lO) C*ii) • (18) 

The density matrix for the system can be written with composite indices, in 
the form 

Pij,ki = —CijCli . (19) 

It has rank one because the system is in a pure state. Now suppose that we 
are performing experiments on one of the subsystems only. Then the relevant 
density matrix is the partially traced density matrix pA = Tr^p, 

n 

ptk = T.P^J^kj . (20) 
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The rank of this matrix may well be greater than one. There are two extreme 
cases. The global state of the system may be a product state, 

n n 

\^>) = \A)\B) = Y.Y.ia,mh,\3)) ^ a,=a.b,, (21) 

i=0 j=0 

so the matrix C is an outer (dyadic) product of two vectors a and b. In this 
case the partially traced density matrix and the matrix Cij both have rank 
one and the individual subsystems are in pure states of their own. A global 
state of this kind is said to be un-entangled or separable. At the opposite 
end of the spectrum it happens that 

1 " 

This means that we know nothing at all about the state of the subsystems, 
even though the global state is precisely known. A global state of this kind 
is said to be maximally entangled. In between these cases are cases where 
the von Neumann entropy of takes some intermediate value. They are 
also entangled. It is known that an arbitrary state in the Hilbert space 
® 'H^ can be brought to the form 

n 

\m) = Y,cm) (23) 

by means of unitary transformations belonging to the subgroup U{N) x U{N), 
that is by "local unitary" transformations acting independently on the two 
subsystems. The coefficients are square roots of the eigenvalues of the density 
matrix obtained when one subsystem is traced out, and hence they obey 

n 
i=0 

This is known as the Schmidt decomposition and is explained in many places 
[|r^|]rT|. Once the Schmidt coefficients q have been ordered (say by their 
sizes) they cannot be changed by local unitaries, so that they can be used 
to label the orbits of U{N) x U{N) in CP'^. This ends our brief review of 
entanglement. 

Now we are going to take a look (literally!) on the sub manifolds that we 
mentioned. Consider first the separable states, as described in eq. (0). If 
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we keep \A) fixed and vary \B) we sweep out a CP" that lies, in its entirety, 
in the submanifold. Reversing the roles of 1^4) and |-B) we sweep out another 
CP", and every state in the submanifold can be reached by a combination 
of these operations. This means that the submanifold of separable states 
is the Cartesian product CP" x CP". For N = 2 (that is n = 1) we get 
the four real dimensional submanifold CP^ x CP^ sitting inside CP^, and 
the question is what it looks like in our picture. We want it to appear as a 
surface in the 3-torus, sitting over a surface in the octant. It is not a priori 
clear that this can be arranged. It depends on the arrangement of the octant 
picture. As a matter of fact it works nicely if the corners of the octant are 
made to represent separable states (and it does not work nicely if the corners 
represent the Bell basis). Explicitly, the state is separable if the rank of Cij 
is unity. Using the homogeneous coordinates introduced already in eq. (|^) 
this means that 

In our coordinates this corresponds to the two real equations 

rioria — nin2 = (26) 

1^1 + 1^2-1^3 = . (27) 

The equation does separate into two equations, one independent of the phases 
and the other involving only them. This is a picture that can be drawn. 

The surface in the 3-torus is, in itself, a flat 2-torus. The surface in the 
octant is easy to draw too because it has an interesting structure. In eq. 
(pTD, keep the state of the second subsystem fixed. Say bo/bi = ke"^, where 
A; is a real number. This implies that 

b 

— = no = kni (28) 

^ = ^ n2 = kns. (29) 

As we vary the state of the first subsystem we sweep out a curve in the 
octant that is in fact a geodesic (the intersection between the 3-sphere and 
two hyperplanes through the origin in the embedding space). In our picture 
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this is a straight hne. In this way we see that the surface in the picture is 
a ruled surface swept out by two famihes of straight hues. If the octant of 
the 3-sphere (that is the interior of our gnomonic tetrahedron) had been flat 
this would have been an intrinsically curved surface, as a glance at Figure 6 
shows. But the true intrinsic geometry of the surface we see in the picture is 
flat. To see this it is convenient to coordinatize the octant by means of Euler 
angles; 



/ no \ 


/ 


ni 








\ ris ) 


V 



Msinl \ 



sm sm 2 
sin ^ cos f 
cos sin I 
cos ^ cos f 



(30) 



A straightforward calculation now shows that the separability condition (^) 
is equivalent to = 0. The coordinate r varies with the state of the first 
subsystem, and the coordinate 6 with the other. The intrinsic metric of the 
surface that we see in the octant is 



idT^ + d9^) . 



(31) 



Hence it is an intrinsically flat surface embedded in a curved space; in the 
language of toric geometry the "manifold with corners" of CP^ x CP^ is a 
flat square. 

This little calculation has an interesting interpretation that we mention 
for the benefit of those readers who are familiar with the Hopf fibration of 
the 3-sphere |]l|. The 3-sphere can be filled with a congruence of nowhere 
vanishing geodesies ("Villarceau circles") that twist around each other but 
never meet. The Euler angles are adapted to this congruence in such a 
way that r runs along the geodesies and the latter are labelled by 6 and 0. 
Thus we see that our surface is made up of a one parameter family of " Hopf 
fibres" ; actually there is another Hopf fibration with the opposite twist and 
this explains why the surface is ruled by two families of intersecting geodesies. 

We are now ready to look at the space of separable states in CP^. It is 
given in Figure 6. 



We now turn to the maximally entangled states. According to eq. (22) a 
state is maximally entangled if and only if the matrix Cij is unitary. Since 
an overall factor of this matrix is irrelevant for the state we can use it to 
adjust the determinant of the matrix to equal one. The state now determines 
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Figure 6: The separable states, or the submanifold CP^ x CP^ in CP^. It 

appears as a surface ruled by geodesies in the octant, and it must be remembered 
that there is a 2-torus sitting in the 3-torus over each of its interior points. The 
intrinsic geometry of the " manifold with corners" , that is the ruled surface that 
we see in the octant, is actually flat. Two different perspectives of the octant are 
shown. The 3-torus is shown schematically as a cube. 
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the matrix up to multiplication with an A^th root of unity. We arrive at 
the conclusion that the space of maximally entangled states form the group 
manifold SU{N)/Z^ = U{N)/U{1) [g]. The reason why this space 



turns up is that, like the separable states, the maximally entangled ones 
forms an orbit of the group of local unitary transformations. The group 
manifold SU{N)/Z^ and the product space CP" x CP" then spring to 
mind as the two most obvious candidates. 

When = 2 we have SU{2)/7P' = S0{3), and this happens to be the real 
projective space RP'^. To see what it looks like in the picture we parametrize 
the unitary matrix Cij as 



a, = j Z^ = («, -P*, a*) . (32) 

In our coordinates this yields three real equations, namely 

no = ns ni = 712 (33) 

Z/l + Z/2 - Z/3 = TT . (34) 

In the octant this is a single geodesic connecting the entangled edges and 
passing through the center of the tetrahedron. In the 3-torus it is a surface 
representing a flat 2-torus. So it appears as a one parameter family of 2-tori, 
in fact an RP^, as it should. Note that the geodesic in the octant crosses 
the separable surface where the 3-torus achieves its maximum size; this is 
how it manages to keep its distance from the entangled states (namely, every 
maximally entangled state is at the distance ir/A from the separable surface). 
This is shown in Figure 7, where we also show the location of the maximally 
entangled Bell basis 

|^±) = |0)|1)±|1)|0) |</>±) = |0)|0)±|1)|1). (35) 

Finally, let us illustrate the collapse of the wave function. If we perform 
a measurement on one of the subsystems the state of the composite system 
will "jump" in the general direction of the separable surface. If the measure- 
ment is a von Neumann measurement we end up on the surface of separable 
states. The question is, where? The most likely possibility is on that point 
on the separable surface that is closest to the state we started out from. In 
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Figure 7: The maximally entangled states form an RP'^'s worth of points, all 
at distance 7r/4 from the closest separable states. In the octant we see a single 
straight line; the location of the standard Bell basis is also shown. There is a 
2-torus in the 3-torus (again shown schematically as a cube) over each interior 
point. 



(0.1.0.0) 




iO.O.I.O) 




(1.0.0.0) 



Figure 8: If the global system is in the Bell state l^*") and if we perform a 
complete measurement on one of the subsystems we end up on a sphere that lies 
on the separable surface. It is a one parameter family of circles as it should be. 
We show one such circle lying on the separable surface in the 3-torus on the right 
hand side. 
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the generic case this point is unique. (To compute it, first use local unitaries 
to bring the given state to the Schmidt form. Unitary transformations are 
isometries and do not change distances. It is now easy to see that the closest 
separable state is in fact the nearest corner of the Schmidt simplex. Gener- 
ically this is unique.) Maximally entangled states are exceptional in that 
there is an entire 2-sphere's worth of points on the separable surface at equal 
distance from the given state. This 2-sphere is not a complex projective line 
because its radius is 1/ -\/2 whereas a complex projective line is a sphere with 
radius 1/2. We are equally likely to land anywhere on this 2-sphere. This is 
basically the statement that if the global state is completely entangled then 
we do not know anything about the state of a subsystem. So Figure 8 picture 
depicts that 2-sphere of separable states that is closest (at distance 7r/4) to 
the Bell state \^~). It consists of all states of the form \+)n\—)n for some 
choice of direction n. 



4. Submanifolds of fixed entanglement. 

We can go on to ask about states with some intermediate degree of entangle- 
ment. In the particular case when N = 2 the remainder of CP^^~^ is foliated 
by a one-parameter family of five dimensional orbits of local unitaries. We 
will not draw such a submanifold here but it can actually be done, and the re- 
sult is much nicer than one might have thought. A pure state of intermediate 
entanglement has the Schmidt decomposition 

1^') = cos(t| -F) -F-sinal ), < a < 7r/4 . (36) 

The entanglement grows monotonically with the Schmidt angle a. The par- 
tially traced density matrix (that arises when we trace over one of the sub- 
systems) then has two unequal and non-zero eigenvalues Ai = cos^ a and 
A2 = singer. A minor calculation verifies that the partially traced density 
matrix has these eigenvalues if and only if 

nlnj + nlnl - cos^ a sin^ a 

cos (i/3 - i/1 - U2) = . (37) 

2nonin2n3 

When this equation has any solutions at all, it describes a 2-torus in the 
3-torus, but this time its position within the 3-torus depends on where we 



20 



are in the octant. There are three dimensions left to account for, and the 
pleasant surprise is that this appears as a volume that only partly fills the 
octant. Its boundaries are obtained by setting the right hand side of the 
preceding equation equal to ±1, that is by 

{nQU^ — 711^2)^ < cos^ a sin^ a < (rio^a + ?t-i?t-2)^ • (38) 

When a is small it lies close to the separable surface, while for a close to 7r/4 
what we see in the octant is a kind of three dimensional tube surrounding the 
maximally entangled line. We leave its precise appearance as an exercise for 
the reader. Note that the case = 2 is exceptional; when N > 2 the orbits 
of local unitary transformations have codimension larger than one, there is 
no obviously canonical measure of pure state entanglement, and a much more 
intricate picture emerges. For illustrations of the = 3 case consult 

The situation is qualitatively similar to that of the orbits of SU{2) in 
CP^ m , also in the sense that the octant picture of the orbit looks nice but 
is poorly suited to do calculations. Once it is understood what a f/(2) x 
U{2) orbit looks like we can watch it grow and shrink as we increase the 
entanglement. But to understand its intrinsic geometry it is better to proceed 
as follows: Choose a point in the orbit given in Schmidt form by 

where Cij is the matrix introduced in eq. (|l^) and a runs from to 7r/4; cr is 
the Schmidt angle and increases monotonically with the entanglement. An 
arbitrary point in the orbit can be reached from the given one by means of 
unitary transformations of the factor Hilbert spaces. Actually SU{2) trans- 
formations are enough, and these can be parametrized with Euler angles. 
That is to say that an arbitary point in the orbit can be given by a matrix 
Cij defined by 

(J — ^-i'i'iLz ^iOiLy ^-iTiL^ ^-iT2Lz ^-i92Ly ^-i(j)2L^ (4:0) 

where Ly and are the angular momentum operators in the standard repre- 
sentation (and one of the group elements appears transposed in the formula). 
In the calculation one sees that only r = ri + r2 matters. After a still fairly 
elaborate calculation we obtain the Fubini-Study metric in the form 
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ds^ = da^ + df , (41) 

where dt^ is the intrinsic metric on the five dimensional orbit labelled by 
< cr < 7r/4. Explicitly 



dt^ = ^{d9l + sin^ 9id4>l + del + sin^ ^2^/02 + 



(42) 



+2 sin 2o"(cos r sin 6i sin 92d(f)id(f)2 — d9id92 — 
— sin r sin 92d9id(f)2 — sin r sin 9id92d(f)i) 

+ sin^ 2a{dT + cos ^^irf^i + cos ^^2*^02)^) 

where < (j)i,(j)2,T < 2tt, < 6^1, 6'2 < vr. When a = this reduces to 
the metric of x as it should, while the coordinates misbehave when 
a = tt/A. The square root of the determinant of this metric is 

y/g = — cos^ 2cr sin 2a sin 9i sin 6*2 . (43) 

(Actually the clever way to compute this is to go via the symplectic form 
presented in section 6. The coordinates used here are well adapted for this 
task.) The volume of a given orbit can now be computed and is found to be 

vol(cr) = TT^ cos^ 2a sin 2a . (44) 

Dividing by the volume vr^/G of CP^ we obtain a probability distribution 
P{a) for the Schmidt angle a. Note that the unitarily invariant distribution 
over the set of pure states used here induces the unitary distribution inside 



the Bloch ball for the mixed states pa = T^''^b\'^){'4^\- 111 • 

Another coordinatization of the orbit works better when the entanglement 
is close to maximal. Above we used two arbitrary unitary matrices ui and 
U2 to write 

C = uiC\2 ■ (45) 
If we introduce us = uiU2 this becomes 
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C = uiC\^^U3 . (46) 

Now one Euler angle in ui is irrelevant; moreover we see directly that in 
the maximally entangled case — when becomes diagonal — the orbit col- 
lapses to the group manifold of SU{2)/Z'^ = RPl One can check that the 
embedding is isometric. 

Since we have a foliation of CP^ with five dimensional hypersurfaces — 
except for the exceptional subsets where the entanglement either vanishes or 
is maximal — it is interesting to look at the second fundamental form of these 
hypersurfaces. One finds that the trace of the extrinsic curvature tensor is 

K = (cos^ 2a -2 sin^ 2a) . (47) 

cos 2a sin 2a 

A foliation with the property that K is constant on each hypersurface is 
called a constant mean curvature foliation. We observe that K = when 
tan 2(7 = I/V2, that is when the volume of the orbit is maximal. 

Let us emphasize again that this state of affairs is special to two qubit 
entanglement; in higher dimensions a much more intricate story emerges. See 
ref. []TB|] for the dimensions of the local orbits that appear in the N x N case. 



5. An Aside on the Statistical Geometry of Simplices. 

Due to obvious limitations we cannot literally draw our picture when the 
dimension of the Hilbert space exceeds four. There is one particular aspect 
of entanglement that we can draw up to dimension 4x4, however. This is 
the Schmidt simplex, that is to say all states that have the form given in 
eq. (|23|). As we mentioned a central fact is that given any state there is a 
local unitary transformation that brings it into the Schmidt simplex. It is a 
simplex essentially because of eq. (p^. It is clear that the Schmidt simplex 
is an — 1 dimensional (hyper-)face of the hyperoctant that forms a part of 
our picture. It is easy to check that its intrinsic geometry — induced by the 
Fubini-Study metric — is round, so in itself it forms the (hyper-)octant of a 
round sphere. Using gnomonic coordinates we draw it as a flat simplex. 

The Schmidt simplex should not be confused with another simplex having 
the same pure states sitting in its corners, namely the statistical simplex 
whose points represent density matrices in the convex cover of these states. 



23 



The latter is an — 1 dimensional subset of the iV^ — 1 dimensional space of 
density matrices, consisting of density matrices that can be simultaneously 
diagonalized. Usually the probabilities Pi are used as barycentric coordinates 
on this simplex, which then consists of all density matrices of the form 

n n 

P = Y.Pi\^)(^\ ^ Ep* = i- (48) 

1=0 i=0 

There is an obvious flat metric on this simplex that has the property that the 
density matrices obtained by taking the mixture of two density matrices in 
the set appears as a straight line connecting two extreme points representing 
the original pair of density matrices. There is also a round metric on this 
simplex that is induced by the Bures metric on the set of all density matrices 
[0; it captures the statistical geometry of density matrices |I8| . 

A possible confusion now arises because we have one round and two flat 
metrics on the same simplex, the latter two being the metric that makes 
statistical mixtures appear as straight lines and the metric that naturally 
exists on the gnomonic coordinate plane. These are not the same. If we 
draw the round simplex using gnomonic coordinates the statistical mixtures 
will be represented by curves, as in Figure 9. This must be kept in mind 
if we ask (say) how far away the corners are from the interior: The usual 
statement P| that the corners are further away than they seem refers to the 
coordinates pf, the significance of these coordinates is that they manifest the 
convexity properties of the simplex. 

Note also that we have two different round simplices with pure states in 
its corners, namely the Schmidt simplex that consists of pure states, and the 
statistical simplex that consists of density matrices. 

6. An aside on symplectic geometry. 

The manifold CP" does not enjoy a metric geometry only. It also has a 
symplectic structure closely allied to it; in physicist's language it is a phase 
space equipped with a Poisson bracket in a natural way. In more mathemat- 
ical language there is a symplectic (closed, non-degenerate) 2-form around. 
Since we are praising the virtues of a special coordinate system here it seems 
natural to mention that the symplectic 2-form Q also takes a simple form in 
this coordinate system. In effect 
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Figure 9: The statistical simplex of diagonal density matrices of size = 3, 
drawn with the eigenvalues pi as barycentric coordinates (left), and the same 
simplex drawn as a round simplex using gnomonic coordinates x, y (right). Corners 
represent pure states while the centers of the triangles represent the maximally 
mixed state. The round (statistical) distance from a given point to the nearest 
corner appears too short in the former picture and too long in the latter. 



Z ■ ZdZ ■ AdZ - Z ■ dZ AdZ ■ Z 

n = 2i = = = 

Z-ZZ-Z 

(49) 

= ^{nxdnx A dui + n2dn2 A dv2 + n^dn^ A du^) , 

where the last line is for n = 3. Hence we can think of as being "canon- 
ically conjugate" to the phase uf, in effect these are action-angle variables. 
There is a nice interplay between the symplectic geometry and the geometry 
of entanglement. In particular a straightforward calculation verifies that the 
space of maximally entangled states is a Lagrangian submanifold. A sub- 
manifold of dimension Z) in a symplectic space of dimension 2D is said to 
be Lagrangian if and only if the symplectic form induced on the subman- 
ifold by the embedding vanishes; mathematicians know that SU{N)/Z^ is 
a Lagrangian submanifold of CP^^~^ but no complete classification of La- 



grangian submanifolds is available |19 . 
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Let us now specialize to CP'^, and give the symplectic form in the coor- 
dinates introduced in eqs. (|39D -(^0D: 



Q = 2 sin 2a{dcr A dr + cos Oidcr A d(f)i + cos 62d(J A d(f)2) + 

(50) 

+ cos 2o"(sin OidOi A d(j)i + sin 92d62 A d(j)2) ■ 

On the five dimensional orbits of local unitaries the first line goes away since 
a is constant. The symplectic form is then degenerate and goes smoothly 
over to the symplectic form on x as a goes to zero. In the language 



often used by physicists interested in constrained systems ||2^ the equation 
( p?!) that defines the orbit is a first class constraint and the coordinate r runs 
along the gauge orbits. 

We observe that the volume element discussed in section 5 is easy to 
derive using 

V^d'x = A i\n) A i\n) . (51) 

This alternative way of computing the volume element is always open on 
Kahler manifolds such as CP". 



7. Envoi. 

We hope that our results serve to bring home the fact that CP"" has an 
interesting geometry, that this geometry can be visualized without too much 
effort, and that the resulting picture has something to tell us about the 
physics of entanglement — at the very least, that it serves to illustrate it. 
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